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A COMPARATIVE ANALYSIS OF THE PERFORMANCE OF 
LONG-RANGE HYPERVELOCITY VEHICLES 1 

By Alfred J. Eggers, Jr., H. Julian Allen, and Stanford E. Neice 


SUMMARY 

Long-range hyj)er velocity vehicle s are studied in terms of their 
motion in powered flight , and their motion and aerodynamic 
heating in unpowered flight . Powered flight is analyzed for 
an idealized propulsion system which approximates rocket 
motors. Unpowered flight is characterized by a return to earth 
along a ballistic , skip , or glide trajectory. Only those trajectories 
are treated which yield the maximum range for a given velocity 
at the end of powered flight. Aerodynamic heating is treated 
in a manner similar to that employed previously by the senior 
authors in studying ballistic missiles (NACA Rep. 1381), 
with the exception that radiant as well as convective heat transfer 
is considered in connection with glide and skip vehicles. 

The ballistic vehicle is found to be the least efficient of the 
several types studied in the sense that it generally requires the 
highest velocity at the end of powered flight in order to attain a 
given range. This disadvantage may be offset, however, by 
reducing convective heat transfer to the re-entry body through 
the artifice of increasing pressure drag in relation to f riction 
drag — that is, by using a blunt body. Thus the kinetic energy 
required by the vehicle at the end of powered flight may be 
reduced by minimizing the mass of coolant material involved. 

The glide vehicle developing lift -drag ratios in the neighbor- 
hood of and greater than 4 is far superior to the ballistic vehicle 
in ability to convert velocity into range. It has the disadvantage 
of having far more heat connected to it; however, it has the 
compensating advantage that this heat can in the main be 
radiated back to the atmosphere. Consequently, the mass of 
coolant material may be kept relatively low. 

The skip vehicle developing lift-drag ratios from about 1 to 4 
is found to be superior to comparable ballistic and glide vehicles 
in converting velocity into range. At lift-drag ratios below 1 it 
is found to be about equal to comparable ballistic vehicles while 
at lift- drag ratios above 4 it is about equal to comparable glide 
vehicles. The skip vehicle experiences extremely large loads, 
however, and it encounters most severe aerodynamic heating. 

^4$ a final performance consideration, it is shown that on the 
basis of equal ratios of mass at take-off to mass at the end of 
powered flight , the hypervelocity vehicle compares favorably 
with the supersonic airplane for ranges in the neighborhood of 
and greater than one half the circumference of the earth. In the 
light of this and previous findings, it is concluded that the 
ballistic and glide vehicles have, in addition to the advantages 


usually ascribed to great speed, the attractive possibility of pro- 
viding relatively efficient long-range flight . 

Design aspects of manned hypervelocity vehicles are touched 
on briefly. It is indicated that if such a vehicle is to develop 
relatively high lift-drag ratios , the wing and tail surfaces should 
have highly swept, rounded leading edges in order to alleviate 
the local heating problem with minimum drag penalty. The 
nose of the body should also be rounded somewhat to reduce 
local heating rates in' this region. If a manned vehicle is de- 
signed for global range flight, the large majority of lift is ob- 
tained from centrifugal force, and aerodynamic lift-drag ratio 
becomes of secondary importance while aerodynamic heating 
becomes of. primary importance. In this case a glide vehicle 
which enters the atmosphere at high angles of attack, and hence 
high lift, becomes especially attractive with a more or less 
rounded bottom to minimize heating over the entire lower surface. 
The blunt ballistic vehicle is characterized by especially low 
heating, and it too may be a practical manned vehicle for ranges 
in excess of semiglobal if great care is taken in supporting the 
occupant to withstand the order of 10 g’s maximum deceleration 
encountered during atmospheric entry. 

INTRODUCTION 

It is generally recognized that hypervelocity vehicles are 
especially suited for military application because of the great 
difficulty of defending against them. It is also possible 
that for long-range operation, hypervelocity vehicles may 
not be overly extravagant in cost. A satellite vehicle, for 
example, can attain arbitrarily longTangc with a finite speed 
and hence finite energy input. E. Sanger was among the 
first to recognize this favorable connection between speed 
and range (ref. 1) and was, with Bredt, perhaps the first to 
exploit the speed factor in designing a long-range bomber 
(ref. 2). This design envisioned a rocket-boost vehicle 
attaining hypervelocities at burnout and returning to earth 
along a combined skip-glide trajectory. Considerable at- 
tention was given to the propulsion and motion analysis; 
however, little attention was given to what is now con- 
sidered to be a principal problem associated with any type 
of hypersonic aircraft, namely that of aerodynamic heating. 
In addition, the category of expendable vehicles, perhaps 
best characterized by the ballistic missile, was not treated. 

Since the work of Sanger and Bredt there have been, of 
course, many treatments of long-range hypervelocity velii- 
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cles in which the propulsion, motion, and heating problems 
have been studied in considerable detail. However, these 
analyses have been devoted in the main to particular designs 
and are not intended to reveal, for example, the relative ad- 
vantages and disadvantages of ballistic-, skip-, and glide- 
type vehicles. Furthermore, it appears that the extent to 
which these vehicles can compete on a simple efficiency basis 
with lower speed aircraft of either the expendable or non- 
expendable type has not been well established. 

It has therefore been undertaken in the present report to 
make a comparative analysis of the performance of hyper- 
velocity vehicles having ballistic, skip, and glide trajectories. 
An idealized propulsion system, whose performance approxi- 
mates that of rocket motors, is assumed. The motion 
analysis is simplified by treating, for the most part, only 
optimum trajectories yielding the maximum range for 
given initial kinetic energy per unit mass in the unpowered 
portion of flight. Aerodynamic heating is treated in a man- 
ner analogous to that employed by the senior authors in 
studying ballistic missiles (ref. 3) with the exception that 
radiant heat transfer, as well as convective heat transfer, is 
considered in the treatment of glide and skip vehicles. The 
efficiencies of these vehicles are compared with supersonic 
aircraft with typical air-breathing power plants. 

NOTATION 

A reference area for lift and drag evaluation, sq ft 

c specific heat of vehicle material, ft-lb/slug °R 

C D drag coefficient 

C L lift coefficient 

C F skin-friction coefficient 

C/ equivalent skin-friction coefficient (see eq. (40)) 

C v specific heat of air at constant pressure, ft-lb/slug 

°R 

C v specific heat of air at constant volume, ft-lb/slug 

°R 

D drag, lb 

e Naperian logarithm base 

E performance efficiency factor (see eq. (85)) 

/ general functional designation 

F a ,F b functions of A J, (see eqs. (74) and (80)) 

G ratio of maximum deceleration to gravity 

acceleration (32.2 ft/sec 2 ) 

g acceleration due to force of gravity, ft/sec 2 

h convective heat-transfer coefficient, ft-lb/ft 2 sec 

°R 

H convective heat transferred per unit area (unless 

otherwise designated), ft-lb/ft 2 
I specific impulse, sec 

J range parameter for glide vehicle (see eq. (68)) 

k Stefan-Boltzmann constant for black body 

radiation (3.7X10 -10 ft-lb/ft 2 sec °R 4 ) 

K constant in stagnation point heat-transfer equa- 

tion, slug 1/2 /ft (see eq. (44)) 

L lift, lb 

m mass, slugs 

M Mach number 

Q convective heat transferred (unless otherwise 

designated), ft-lb 


/* distance from center of the earth, ft 

r e radius of curvature of flight path, ft 

r 0 radius of earth, ft 

R range, ft 

$ distance along flight path, ft 

5 surface area, sq ft 

t time, sec 

T temperature (ambient air temperature unless 

otherwise specified), °R 

V velocity, ft/sec 

V ratio of velocity to satellite velocity 

V s velocity of satellite at earth’s surface (taken as 

25,930 ft/sec) 

W weight, lb 

y vertical distance from surface of earth, ft 

a angle of attack, radians unless otherwise speci- 

fied 

j3 constantin density-altitude relation, (22,000 ft” 1 ; 

see eq. (15)) 

7 ratio of specific heats, C v /( \ 

6 semivertex angle of cones, radians unless other- 

wise specified 
A increment 

rj lift-drag efficiency factor, (see eq. (B27)) 

6 angle of flight path to horizontal, radians unless 

otherwise specified 

A leading edge sweep angle, deg 

p air density, slugs/cu ft (p 0 = 0.0034) 

<j nose or leading-edge radius of body or wing, ft 

partial range, radians 
3> total range, radians 

^ remaining range (4>— <p), radians 


Subscripts 

conditions at zero angle of attack 
1, 2, 3, . . . conditions at end of particular rocket stages 
a conditions at point of maximum average heat- 

transfer rate 
av average values 

b conditions at point of maximum local heat- 

transfer rate 
C convection 

e effective values 

en conditions at entrance to earth’s atmosphere 

ex conditions at exit from earth’s atmosphere 

f conditions at end of powered flight 

i initial conditions 

l local conditions 

n ballistic phases of skip vehicles 

N total number of rocket stages 

p pressure effects 

P pay load 

r recovery conditions 

R radiation 

« stagnation conditions 

T total values 

w wall conditions 
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ANALYSIS 


GENERAL CONSIDERATIONS 


In the following analysis of long-range hypervelocity 
vehicles, only flight in planes containing the great circle 
arc between take-off and landing is considered. The flight 
is thought of in two phases: (a) the powered phase in which 
sufficient kinetic energy, as well as control, is imparted to the 
vehicle to bring it to a prescribed velocity, orientation, and 
position in space; and (b) the unpowered phase, in which the 
vehicle travels to its destination under the influence of 
gravity and aerodynamic forces. 

The analyses of motion and aerodynamic heating during 
unpowered flight will, of necessity, differ widely for the 
several types of vehicles under consideration. On the other 
hand, motion in the powered phase is conveniently treated 
by a method common to all vehicles. The study of powered 
flight and its relation to range is therefore taken as a starting 
point in the analysis. 


POWERED FLIGHT AND THE BREGUET RANGE EQUATION 

In this part of the study, the following simplifying as- 
sumptions are made: (a) aerodynamic heating can be 
neglected on the premise that high flight speeds are not 
attained until the vehicle is in the rarefied upper atmosphere; 2 

(b) sufficient stability and control is available to provide 
proper orientation and positioning of the vehicle in space; 

(c) ^the distance traveled while under power is negligible by 
comparison to the overall range; and finally, (d) the thrust 
is very large compared to the retarding aerodynamic and 
gravity forces. In terms of present-day power plants, the 
last assumption is tantamount to assuming a rocket drive 
for the vehicle. 

The velocity at burnout of the. first stage of a multistage 
rocket (or the final velocity of a single-stage rocket) can then 
be expressed as (see, e. g., ref. 4) : 



where the initial velocity is taken as zero. In this expression, 
ra t and m fl represent the mass of the vehicle at the beginning 
and ending of first-stage flight, and V fl = Y fl /V s where V s = 
^gr 0 = 25,930 feet per second is the satellite velocity at the 
surface of the earth. The coefficient g is the acceleration due 
to gravity and is, along with the specific impulse 7, con- 
sidered constant in this phase of the analysis. The final 
velocity of the vehicle at the end of the N stages of powered 
flight can be expressed as 




where the initial mass of any given stage differs from the 
final mass of the previous stage by the amount of structure, 
etc., jettisoned. 

Now let us define an equivalent single-stage rocket having 
the same initial and final mass as the A r -stage rocket and the 


same initial and final velocity. There is, then, an effective 
specific impulse defined by 




whereby equation (2) can be written as 

v '=& in (s-;) <4) 


The effective specific impulse I e is always somewhat less 
than the actual specific impulse, but for an efficient design 
they are not too different. Throughout the remainder of 
the analysis the effective impulse I e will be used. 

Equation (4) might be termed the “ideal power plant” 
equation for accelerated flight because, when considered in 
combination with the assumptions underlying its develop- 
ment, attention is naturally focused on the salient factors 
leading to maximum increase in velocity for given expendi- 
ture of propellant. Thus the thrust acts only in over- 
coming inertia forces, and the increase in vehicle velocity 
is directly proportional to the exhaust velocity (gl) of the 
propellant. 

Now we recognize that an essential feature of the hyper- 
velocity vehicles under study here is that they use their 
velocity (or kinetic energy per unit mass) to obtain range. 
For this reason, equation (4) also constitutes a basic per- 
formance equation for these vehicles because it provides 
a connecting link between range requirements and power- 
plant. requirements. 

In addition to comparing various types of liypervelocity 
vehicles, our attention will also be focused upon comparison 
of these vehicles with lower speed, more conventional types 
of aircraft. For this purpose it is useful to develop an 
alternate form of equation (4). We observe that the 
kinetic energy imparted to the vehicle is 


hrrifV / 


This energy is equated to an effective work done, defined as 
the product of the range traveled and^ a constant retarding 
force. (Note that the useful kinetic energy at the end of 
powered flight is zero.) This force is termed the “effective 
drag” D e . Thus 

D e R=% m f V f 2 (5) 


where R is flight range measured along the surface of the 
earth. Similarly, we may define an “effective lift” L e , 
equal to the final weight of the vehicle 

L e =W f =m f g 


from which it follows that equation (5) may be written as 



( 6 ) 


where ( L/D) e is termed the “effective lift-drag ratio.” 
Combining equations (4) and (6) , we obtain 





* This assumption is in the main permissible. A possible exception occurs, however, with 
the glide vehicle for which heat-transfer rates near the end of powered flight can be comparable 
t o those experienced in unpowered gliding flight 


( 7 ) 
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where 



( 8 ) 


and represents an “effective” flight velocity of the vehicle. 
Equation (7) will prove useful in comparing hypersonic 
vehicles with conventional aircraft because of its analogy 
to the Breguet range equation, 

*=b /f 1 "(S;) <*> 


It will also prove useful to have equation (7) in the dimen- 
sionless form obtained by dividing through with r 0} the 
radius of the earth. In this case we have 



( 10 ) 


where $ is the range in radians of arc traversed along the 
surface of the earth. 


MOTION IN UNPOWERED FLIGHT 

Ballistic trajectory. — In studying the motion of long-range 
vehicles in this trajectory, advantage is taken of the fact 
that the traverse through the earth’s atmosphere generally 
forms only a small part of the total trajectory. Therefore, 
the deflection and deceleration encountered in the re-entry 
phase (discussed in detail in ref. 3) are neglected in the 
computation of the total range and rotation of the earth is 
neglected in this and all other phases of the analysis. With 
the added simplification that the contribution to range of 
the powered phase of flight is negligible, the ballistic tra- 
jectory becomes one of Kepler’s planetary ellipses, the 
major axis of which bisects the total angle of arc <t> traveled 
around the earth. For the trajectories of interest here 
(E/<1), the far focus of the ellipse is at the mass center 
of the earth. For purposes of range computation, then, the 
ballistic vehicle leaves and returns to the earth’s surface 
at the same absolute magnitude of velocity and incidence 
(see sketch) . 



The expression for range follows easily from the equation 
of the ellipse (see, e. g., ref. 5) and can be written 


sin 6 f cos d\ 

=U«»'J <>» 

v,‘ 7 

where the angle of incidence 0 r is considered positive. In 
order to determine the optimum trajectory giving maximum 
range for a given velocity V h equation (11) is differentiated 
with respect to 0 f and equated to 0, yielding 


$=-=2 tan" 1 
ro I 


1 7 2= i 72 = 1 — tan 
V s 

<*>=--=»•— 40, 

■> *0 


( 12 ) 


Equations (11) and (12) have been employed to determine 
velocity as a function of incidence for various values of 
range and the results are presented in figure 1. The 'mini- 
mum velocity line” of figure 1 corresponds to the optimum 
trajectories (eqs. (12)). 

The effective lift-drag ratios can easily be calculated for 
optimum ballistic vehicles using equation (6) in combination 
with the information of figure 1. The required values of 
(L/D) e as a function of range are presented in figure 2. 

Skip trajectory. — This trajectory can be thought of as a 
succession of ballistic trajectories, each connected to the 
next by a ‘ ‘skipping phase” during which the vehicle enters 
the atmosphere, negotiates a turn, and is then ejected from 
the atmosphere. The motion analysis for the ballistic 
missile can, of course, be applied to the ballistic phases of 
the skip trajectory. It remains, then, to analyze the 
skipping phases and to combine this analysis with the bal- 
listic analysis to determine over-all range. 

lo this end, consider a vehicle in the process of executing 
a skip from the atmosphere (see sketch). 


+ 

/ 



The parametric equations of motion in directions perpen- 
dicular and parallel to the flight path s are, respectively, 


mV 2 


pV 2 

C L —n~ A— mg cos 0= 

1 r c 

n P y2 a, ■ a dV 
—C D ~ 2 ~ A+mg sin 0=m 


(13) 


where r e is the local radius of curvature of the flight path, 0 is 
the local inclination to the horizontal (positive downward), 
p is the local air density, and C L and C D are the lift and drag 
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Figure 1 .— Variation of velocity with incidence angle for various 
values of range of ballistic vehicle. 


coefficients, respectively, based on the reference area, A, of 
the aircraft. 

In the turning process, aerodynamic lift must obviously 
predominate over the gravity component, mg cos 0. By anal- 
ogy to the atmospheric re-entry of ballistic missiles (see ref. 
3), aerodynamic drag generally predominates ever the gravity 
component, mg sin 0. Moreover, the integrated contribution 
to velocity of this gravity component during descent in a skip 
is largely balanced by an opposite contribution during ascent. 



Figure 2. — Variation of effective lift-drag ratio with range for optimum 
ballistic vehicle. 



Figure 3. — Trajectory of the first skipping phase for a skip vehicle 
with a lift-drag ratio of 2 and a total range of 3440 nautical miles 
(*= 1 ). 


For these reasons we will idealize the analysis by neglecting 
gravity entirely. This approach is analogous to the classical 
treatment of impact problems in which all forces exclusive of 
impact forces (aerodynamic forces in this case) are neglected 
as being of secondary importance. Gravity is shown to be 
of secondary importance in figure 3 where the trajectory re- 
sults obtainable from equations (13) and (14) are presented 
for the first skipping phase of an L/D= 2, 4>=1 skip missile. 

With gravity terms neglected, equations (13) reduce to 


\c LP V'A=-mV'f s 

1 dV 

~2 ^A=m f 


04) 


where dd/ds= to the accuracy of this analysis. 

r c 

Now we assume an isothermal atmosphere, in which case 


p=p 0 e & 


(15) 


where p 0 and p are constants, and y=(r— r 0 ) is the altitude 
from sea level (see ref. 3 for discussion of accuracy of this as- 
sumption). Noting that dy/ds= —sin 0, we combine the first 
of equations (14) with equation (15) to yield 


e~ Bv dy = sin 0 dd 


06 ) 


This expression can be integrated to give 


ClPqA 

2 pm 


-/9V- 


= COS 0 — COS I 


(17) 


where o is t a as zero at the altitude corresponding to the 
effective “outer reach” of the atmosphere. Equation (17) 
points out an important feature of the skip path; namely, 
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cos 6 is a single-valued function of altitude. Since 6 proceeds 
from positive to negative values, it is evident that 

( 18 ) 

where the subscripts en and ex refer to atmospheric entrance 
and exit conditions, respectively, and the numbers n— 1 and 
n refer to successive ballistic phases of the trajectory. Now 
since 

dV _ydV_\ dV 2 
dt ds 2 ds 


equations (14) may be combined to obtain 

1 dV> = V 2 de 

2 ds LID ds 

which, for constant L/D, can be integrated to yield 


(19) 


V t 


=e lid 


( 20 ) 


en n - 1 


With the aid of equation (18), this expression may be 
written 

( 21 ) 


V ei 


-=e L ' D 


which relates the velocities at the beginning and end of a skip 
to the lift-drag ratio and the entrance angle of the vehicle to 
the earth's atmosphere. From equation (18) it follows 
further that the entrance angle for each skip in the trajectory 
is the same, so that 

0«V=0«» m _i= • • • = 0 / 

and hence equation (21) becomes 





( 22 ) 


We now combine this result of the skip analysis with that 
of the ballistic analysis to obtain the total flight range. 
From equation (11) the range of the nth ballistic segment of 
the trajectory is 


<p n = 2 tan -1 


sjn Of cos 6 f 


(23) 


Consistent with the idealization of the skipping process as an 
impact problem, we neglect the contribution to range of each 
skipping phase so that the total range is simply the sum of 
the ballistic contributions. From equations (22) and (23) 
this range is then 


$=;r=Z>»=2Z) tan" 1 

• o n=l n=l 


sin d f COS dr 


4(w-l)fl f 

z=L- e L,D —cos 2 6 f 

w 


(24) 


which maximizes the range of an aircraft developing a 
particular lift-drag ratio. These skipping angles have been 
obtained with the aid of an IBM CPC, and the corresponding 
values of V f as a function of range for various L/D are 
presented in figure 4. Corresponding values of ( L/D) e have 
been obtained using equation (6) and the results are shown 
in figure 5. 



Figure 4. — Variation of velocity with range for various values of 
lift-drag ratio for skip vehicle. 


Glide trajectory. — The trajectory of the glide vehicle is 
illustrated in the accompanying sketch. As in the previous 
analyses, the distance covered in the powered phase will be 
neglected in the determination of total range. 


L 


D 



^'Earth's surface 


The parametric equations of motion normal and parallel 
to the direction of flight are the relations of equations (13) 
rewritten in the form 


L—mg cos d= 


mV 2 


n . • * dV 

— D+mg sin d=m -jj 


From this expression we see that for any given velocity 
at the end of powered flight there is a definite skipping angle 


( 25 ) 
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Figure 5. — Variation of effective lift-drag ratio with range for various 
values of aerodynamic lift-drag ratio of skip vehicle. 


Under the assumption of small inclination angle 0 to the 
horizontal (thus cos 0~ 1 , sin 0 — 0), constant gravity accelera- 
tion ^i. e., and noting the following relations 


(lV_ y dV_\ dV 2 ' 
dt ds 2 d.s 

l = rf(*-0) 

r c ds 


(26) 


d\// _cos 0_ 1 
ds r ~ r 0 


equations (25) can he written in the forms 


L=—mV 2 7— 

n 1 dV2 . 

- m — + mgd 


mV 2 


r 0 




y 




(27) 


Dividing the first of equations (27) by the second yields the 
following differential equation 


!, \ D d ) + \2D ds 1 ds) r 0 ~° (28) 


L 


But, as is demonstrated in Appendix A, the terms gO and 


V 2 ^ may be neglected so that equation (28) reduces to 


dV 2 


2 v 2 +-^ik=0 


( 29 ) 


Since 


ds r 0 (L/D) 1 LID~ 

V s *=gr 0 

equation (29) can be integrated for constant to give the 
velocity in nondimensional form as 


2(f 


V 2 =l-(l-V/)e L/D 


(30) 


This expression gives velocity as a function of range for what 
Sanger (ref. 2) has termed the equilibrium trajectory — that 
is, the trajectory for which the gravity force is essentially 
balanced by the aerodynamic lift and centrifugal force, or 


It follows from equation (31) that velocity can be expressed 
in the form 


V 2 


1 


CrAVlp 

^ 2 mg 


(32) 


Now it is intuitively obvious that as the maximum range is 
approached, L/W—> 1 and hence V 2 becomes small compared 
to one (see cq. (31)). In this event it follows from equation 
(30) that the maximum range for the glide vehicle is given by 




(33) 


The relation between velocity and range has been deter- 
mined with equation (33) for various values of L/D and the 
results are presented in figure 6. Corresponding values of 
(L/D) e have been obtained using equation (6) and are 
presented in figure 7. 

These considerations complete the motion analysis and 
attention is now turned to the aerodynamic heating of the 
several types of vehicles under consideration. 



Figure G. — Variation of velocity with range for various values of 


lift-drag ratio of glide vehicle. 


483423—59 2 
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Figure 7.— Variation of effective lift-drag ratio with range for various 
values of aerodynamic lift-drag ratio of glide vehicle. 


HEATING IN UNPOWERED FLIGHT 

General considerations. — Three aspects of the aerodynamic 
heating of hypervelocity vehicles will be treated here; namely, 

1. The total heat input 

2. The maximum time rate of average heat input per unit 

area 

3. The maximum time rate of local heat input per unit 

area 

Total heat input is, of course, an important factor in deter- 
mining over-all coolant weight, whether the coolant be solid 
f (e. g., the structure), liquid, or gas, or a combination thereof. 
The maximum time rate of average heat input per unit area 
can determine peak average flow rates in the case of fluid 
coolants and may dictate over-all structural strength in the 
event that thermal stresses predominate. 

Excessive local heating is, of course, a serious problem with 
hypervelocity vehicles. This problem may vary depending 
upon the type of the vehicle. Thus, for the ballistic vehicle, 
an important local “hot spot” is the stagnation region of the 
nose, while for the skip or glide vehicle attention may also be 
focused on the leading edges of planar surfaces used for de- 
veloping lift and obtaining stable and controlled flight. In 
this analysis attention is, for the purpose of simplicity, re- 
stricted to the “hot spot” at the nose. In particular, we 
consider the maximum time rate of local heat input per unit 
area because of its bearing on local coolant flow rates and 
local structural strength. 

It is undertaken to treat only convective heat transfer at 
this stage of the study. As will be demonstrated, radiant 


heat transfer from the surface should not appreciably in- 
fluence convective heat transfer to a vehicle. Therefore, 
alleviating effects of radiation are reserved for attention in 
the discussion of particular vehicles later in the paper. This 
analysis is further simplified by making the assumptions that 

1. Effects of gaseous imperfections may be neglected 

2. Shock-wave boundary-layer interaction may be neg- 

lected 

3. Prandtl number is unity 

4. Reynolds analogy is applicable 

These assumptions are obviously not permissible for an accu- 
rate quantitative study of a specific vehicle. Nevertheless 
they should not invalidate this comparative analysis which is 
only intended to yield information of a general nature regard- 
ing the relative merits and problems of different types of 
vehicle (see ref. 3 for a more complete discussion of these 
assumptions in connection with ballistic vehicles). 

In calculating convective heat transfer to hypervelocity 
vehicles, the theoretical approach taken in reference 3 for 
ballistic vehicles is, up to a point, quite general and can be 
employed here. Thus, on the basis of the foregoing assump- 
tions, it follows that for large Mach numbers, the difference 
between the local recovery temperature and wall temperature 
can be expressed as 


(T r -rj,==^r (34) 

It is clear, however, that the walls of a vehicle should be 
maintained sufficiently cool to insure structural integrity. 
It follows in this case that the recovery temperature, at 
hvpervelocities will be large by comparison to the wall tem- 
perature and equation (34) may be simplified to read 


Tr 


l 


Z! 

2 C p 


(35) 


To the accuracy of this analysis, then, the convective heat 
transfer is independent of wall temperature. Therefore, as 
previously asserted, radiant heat transfer should not appre- 
ciably influence convective heat transfer and the one can be 
studied independently of the other. 

Now, according to Reynolds analogy, the local heat- 
transfer coefficient hi is, for a Prandtl number of unity, given 
by the expression 

hi=\C Fl C VlPl Vi (36) 

where C Ff is the local skin-friction coefficient based on con- 
ditions just outside the boundary layer. With the aid of 
equations (35) and (36) the time rate of local heat transfer 
per unit area, 

^-=hi(T r -TJ t (37) 


can be written as 


(IH = V^_ 

dt 4 C v 


(C Fl C PlPl V,) 


(38) 
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Equation (38) can be integrated over the surface of a body 
to yield the time rate of total heat input as follows 


wherein C Pl is set equal to C p and 

c ’'-iL c '-7V‘ ls 


(39) 


(40) 


The parameter C/ is termed the “equivalent skin-friction 
coefficient” and will be assumed constant at a mean value 
for a particular vehicle. From equation (39) we can obtain 
two alternate forms which will prove useful; namely, the 
altitude rate of total heat input defined bv (note that dy is 
negative for dt positive) 


dQ_ 


1 dQ_pV 2 C/S 


dy F sin 6 dt 4sin0 

and the range rate of total heat input defined as 

dQ = 1 dQ_ p V*C F 'S 

d(r 0 <p) V cos d dt 4 cos 0 


equations (39), (41) or (42), depending upon the particular 
variable used. 

The time rate of average heat input per unit area may be 
obtained from equation (39) as 


dH a r _1 dQ 1 T/ 3 sy t 

dt ~S dt~ 4 pl (f 


(43) 


Consider next the local convective heat transfer in the 
region of the nose. The time rate of local heat input per 
unit area was determined in reference 3 under the assump- 
tions that viscosity coefficient varies as the square root of 
the absolute temperature, and that flow between the bow 
shock wave and the stagnation point is incompressible. In 
this case it was found that 


d -§- K y[i v ‘ 


(44) 


where Z£=6.8X10 -6 . A more detailed study of stagnation 
region flow, including effects of compressibility and dissoci- 
ation of air molecules (ref. 6), shows that the constant, K , 
should have a value more like twice the above value at the 
hypervelocities of interest here. 

With these relations we are now in a position to study 
the heating of the several types of vehicles of interest. 

Ballistic vehicle. — The heating for this case has already 
been analyzed in reference 3. Only the results will be given 
here. 

The ratio of the total heat input to the initial kinetic 
energy was found to be 


Q 


1 C/S / CpPoA \ 

2 C d A \l-£ * msin *r) (45) 


For the “relatively light missile, ” which is of principal interest 
here, 


G D<‘o A 

e 0m sin 6 , <^<^1 

and equation (45) reduces to 

Q 1 CVS 
%m,V r 2 2 C d A 


(46) 


(47) 


The time rate of average heat input per unit area was 
found to be 


ZC DPo A e -p u 


dHg V __ 1 dQ_ CV Pu \ / ^ _py^ 20m sin 0 f 


(48) 


dt S dt 4 
which has the maximum value 


(41) 

(dUV 

\ (dHV 

V 0 Of' TT 3 ■ n 

\ a ~6e CdA m Vfi sln ^ 

(49) 


\ dt , 

Lx ( dt j 


at the altitude 




(42) 


Va= ^ In 

/ 3 CoPoA \ 
\2 fim sin 9 f ) 

(50) 


Equation (49) applies, of course, only if the altitude, y a , is 
above ground level. If the value of y a is negative then the 
maximum average heating rate will, of course, occur at sea 
level. 

The time rate of local heat input per unit area to the stag- 
nation region of the nose was found to be 


IMVt v ‘ 

having a maximum value of 


_ &y ^CpppA ^ 

? ~2 e ~20m sin Q f “ 


Pu 


(51) 


/ dH,\ _/dH a \ _ K //3m sin 6, X7 , 
V dt ) max \ dt ) Vh V 3eC D aA 1 

occurriiig at the altitude 

y„=~ In ( Z pnPpA \ 

^ /3 V/Smsin^/ 


(52) 


(53) 


If the value of y b is negative, then the maximum value 
occurs at ground level. 

Skip vehicle. — With the aid of equation (17), the density at 
any point in a given skipping phase is found to be 


P—p 0 e ft '= 


2 ()m , 

--qTa C0S cos ^ 


(54) 


where it is to be recalled that d en =6f. The corresponding 
velocity for constant L/D is, from equation (19), 


V=V en e L,D 


(55) 
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By substitution of equations (54) and (55) into equation (-39), 
the time rate of total heat input at any point in a skipping 
phase can be expressed as follows: 

3 {Of-O) 

<56> 

Now, recalling that ds/dt=V, the first of equations (14) may 
be combined with equation (17) to yield 


^=— 0V (cos 9 — cos Of) 
at 

and we note further from equation (22) that 


(Vcnh = (Venh 

(r„),_i v, 


-(71-1) 


=e 


LID 


(57) 


(58) 


With the aid of equations (57) and (58), equation (56) can 
be integrated to give the total heat input for a given skipping 
phase, thus we obtain 

-£) (59) 

§ mV, 2 2 CdA \l-e LID J e 

The total heat input for the entire trajectory can be ob- 
tained by summing up the heat inputs for each separate 
skipping phase. Performing this operation yields 


2 

or 


Q _ 

hmV? 


^ _g„ 1 C/S ( _J*t\ f, 

^ A m V, 2 2 CdA VI —e LID )t=‘x 

Q _1 C/S 

\ mV, 2 2 C d A 


co — (n — 1) 


e 


L/D 


(60) 


(61) 


which is identical to the result obtained for the light ballistic 
missile (eq. (47)). This result applies, in fact, to all h\pei 
velocity vehicles which lose the large majority of their kinetic 
energy during atmospheric entry. 

The time rate of average heat input per unit area is 
obtained by dividing equation (56) with the surface area, 
thus yielding 

3(fl,-») 

~f=\ Tjb V ' n (cos 6 ~ C0S 6f)e L ' D (62) 

It can be shown that this expression has a peak value at a 
point in the skip, 9 a , given by 

L/D . . 

(cos ft,— cos 9 , ) =-y- sin 9 a 


or 


ft,= tan' 


L/D 


-sin' 


cos 6 f 




(631 


From, equation (22) it can be concluded that the maximum 
heat-transfer rate will occur in the first skip where T e „= 1 ,; 
consequently, 


3 (»,-»«) 


(dHaA (64) 

V dt La, 6 CdA ' 


The time rate of local heat input per unit area in the 
stagnation region of the nose is obtained by introducing 
equations (54) and (55) into equation (44) with the following 
result: 




3 (0f-0) 
L/D 


(65) 


Equation (65) has a peak value at a point 0 b in a skip given 

by 


x LID . „ 
(cos Q b — cos 0,) = sin 6 b 


or 


0 ft =tarr 


L/D 


-sin' 


cos Or 




( 66 ) 


It is clear in this case also that the heat-transfer rate will 
have its maximum value in the first skipping phase where 
the velocities are highest. Since V r en = 1 f in the first skip, 
equation (65) becomes 

3 (0 f -0 b ) 

3 L/D (A 7 ) 

1 ’ 

Glide vehicle. — From equations (30) and (32), the density 
at a point in the glide trajectory is found to be 


P= 


2 mg (1 — V f 2 )e J 


where 


ClAVs 2 l — 

2(p 


( 68 ) 


J= 


L/D 


By substitution of equations (30) and (68) into equation 
(39), the time rate of total heat input can be expressed as 


^ (1-V,V[1-(1-V,V] M 


(69) 


If equations (30) and (33) are combined with this expression, 
we again obtain 

Q 


1 C/S 


A mV, 2 2 CdA 


(70) 


for the heat transfer to a hypervelocity vehicle during atmos- 
pheric entry. 

Now the time rate of average heat input per unit area of 
a glide vehicle is found by dividing equation (69) with the 
surface area, thus yielding 

d -~=\ ( i-v,v[i-(i-v,v]* (7i) 

It, follows from this expression that the maximum time rate 
of average heat input per unit area is 


((IHa A (dH„\ _ 1 0/ mg Vs 

V dt La, V dt ) Ja 3 N 3 CdA L/D 

at a value J„ given by 

J a =— In |(1 — V, 2 ) 


(72) 


( 73 ) 
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If J a is taken as a reference value, and equations (71) and 
(72) expressed in terms of J a and incremental changes 
AJ=J—J 0 , it can easily be shown that 


(iHgJdt 
(i dHJdt ) 

max 


e* J (3-2e AJ )X=F a (AJ) 


(74) 


The dependence of F a (AJ) on A J is shown in figure 8. 

The velocity at which the maximum average heat input 
rate occurs can be obtained by substituting equation (73) 
into equation (30) yielding 



(75) 


Equations (72), (73), and (74) apply, of course, only when 

F,>(l/V3). 

For cases when V f <(l/^3 ), the maximum time rate of 
average heat input per unit area will occur at the start of 
unpowered flight and is given by 



A J 

Figure 8. — Variations of F a (Aj) and F h {AJ) with AJ. 

DISCUSSION 


/ iHg / dHg A _1 Cr_ rngV s 

\dt ) m ax V dt )j c _ o 2 r D A LjD 


(1-V/)V> 


(76) 


The maximum time rate of local heat input per unit area 
in the stagnation region of the nose is found by first substi- 
tuting equations (30) and (68) into equation (44) to obtain 

w-a- who-Vw (77) 

The maximum time rate is then 



occurring at a value of J b given by 


(78) 


J*=-ln 3(1-17) 


(79) 


PERFORMANCE OF HYPERVELOCITY VEHICLES 

In this study the point of view is taken that the perform- 
ance of long-range hypervelocity vehicles is measured by 
their efficiency of flight. Thus, for example, it is presumed 
that the advantages (military and otherwise) of short time 
of flight accrue equally to all vehicles. 

The efficiency of flight is perhaps best measured by the 
cost of delivering a given pay load a given range — the higher 
the cost, the lower the efficiency. Quite obviously it is far 
beyond the scope of the present paper to actually compute 
this cost. Rather, then, we adopt a more accessible param- 
eter of hypervelocity flight, namely, the initial mass of the 
vehicle, as a measure of cost. In effect, then, the assumption 
is made that the higher the initial mass of a vehicle the 
higher the cost and the lower the efficiency. With these 
thoughts in mind, it is constructive to reconsider the basic 
performance equation (eq. (4)) written in the form 


With J b as a reference, it can easily be shown that 


= m / e Vf,9l<! 


(83) 


where 


dH s /dt 1 

(dH.ldt) mu T 2 * 


(3— e AJ ) =F t (AJ) 


A J — J — J /, 


(80) 


The dependence of F b (Aj) on A J is shown in figure 8. 

With reference to equations (30) and (74) it can be seen 
that the maximum time rate of local heat transfer in the 
stagnation region occurs when 


V= 



(81) 


It is apparent then that equations (77), (78), a nd ( 79) apply 
only when V f >^2/3. For cases where V f <^ 2/3, the max- 
imum time rate of local heat input per unit area will occur 
at the start of unpowered flight and is given by 


(dHs\ / dHs\ 

V dt / mar \ dt ) j q 


= K V s 'V/(\ - V/)'i (82) 


This expression clearly demonstrates the roles played by the 
three factors which influence the initial mass of a vehicle 
required to travel a given range. For one thing there is the 
power plant, and as we would expect, increasing the effective 
specific impulse increases the over-all efficiency of flight in 
the sense that it tends to reduce the initial mass. The 
velocity at burnout influences initial mass by dictating the 
amount of fuel required, and it is not surprising that de- 
creasing the required burnout velocity (e. g., by increasing 
the L/D of a skip or glide vehicle) tends to decrease the 
initial mass. Finally, we see that the initial mass is propor- 
tional to the final mass which consists of the pay load, 
structure (and associated equipment), and coolant. If we 
presume the mass of the pay load to be some fixed quantity, 
then the initial mass will vary in accordance with this mass 
of structure and voolant. 

Now we assume for comparative purposes that the power 
plant for one vehicle is equally as good as the power plant 
for another vehicle — that is to say I e is a more or less fixed 
quantity. In this event it is permissible to restrict our 
attention to two main performance considerations; namely, 
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the prescribed motion as it influences the required burnout 
velocity, and the resulting aerodynamic heating as it influ- 
ences structure and coolant. We therefore proceed to 
discuss the comparative performance of long-range hyper- 
velocity vehicles in terms of these considerations. 

Motion. — The dependence of burnout velocity F/On range 
was determined in the analysis of motion in unpowered 
flight and the results obtained for the several types of hyper- 
velocity vehicles under study were presented in figures 1, 4, 
and 6. Using these results in combination with the basic 
performance equation we have calculated the corresponding 
initial to final mass ratios mi/m f as a function of range. 
For these and subsequent calculations it has been assumed 
that the rocket power plant develops an effective specific 
impulse of 300 seconds. The results of these calculations 
are presented in figure 9 and we observe that, in general, 
the mass ratios are highest for the ballistic vehicle. The 
glide and skip vehicles have comparable and relatively low 
mass ratios at lift-drag ratios in the neighborhood of 4 and 
greater. The skip vehicle is superior, however, to the glide 
vehicle at lift -drag ratios in the neighborhood of 2. From 
considerations of motion alone, then, we conclude that the 
skip vehicle and the glide vehicle developing lift-drag ratios 
greater than 2 are superior efficiency wise, in the sense of this 
report, to the ballistic vehicle. Let us now determine how 
these observations are modified by considerations of aero- 
dynamic heating. 

Aerodynamic heating. — The analysis has revealed one par- 
ticularly salient factor in regard to the heat transferred by 
convection to hypervelocity vehicles that expend the 
majority of their kinetic energy of flight in traveling through 
the earth’s atmosphere. This factor is that the amount of 
kinetic energy which appears in the body in the form of 
heat is proportional to the ratio of friction force to total 
drag force acting on the body (see eqs. (47), (61), and (70)). 
With the possible exception of the relatively heavy ballistic 
vehicle (see ref. 3) all of the hypervelocity vehicles treated 
here do expend the major part of their kinetic energy in 
flight. It is, in fact, only by virtue of this expenditure of 
energy that the skip and glide vehicles achieve long range. 
From the standpoint, then, of reducing the total heat trans- 
ferred by convection, the problem is to determine how the 
ratio of friction force to total drag force can be reduced. 
This matter was discussed in detail in reference 3 in connec- 
tion with ballistic vehicles and it was demonstrated that the 



Figure 9. — Variation of mass ratio with range for various lift-drag 
ratios of hypervelocitv vehicles. 


ratio could be reduced by employing high-pressure-drag 
(i. e., blunt) shapes. It would be most fortunate if this 
avenue of solution were open also to the skip and glide 
vehicles; however, it is readily apparent that such is not the 
case. This conclusion follows simply from the fact that t lie 
skip and glide vehicles must develop reasonably high lift- 
drag ratios to achieve long range. But, as is well known, 
high lift-drag ratios and high pressure drag are incompatible 
aerodynamic properties. Evidently, then, the skip and 
glide vehicles will be relatively slender and they will, by 
comparison to blunt ballistic vehicles, be required to absorb 
large amounts of their kinetic energy of flight in the form 
of heat. On the basis of the calculations of reference 3, it 
does not seem feasible for slender hypervelocity vehicles to 
absorb and retain so much heat (of the order of one-tenth 
the kinetic energy of flight). We are led, therefore, to con- 
sider the possibility of radiating part or all of this heat back 
to the atmosphere. 

Let us first consider radiation heat transfer from the sur- 
face of a glide vehicle. For purposes of simplicity we pre- 
sume a vehicle conical in shape. The base diameter is taken 
as 3 feet and the weight as 5,000 pounds. W'e consider two 
slender cones which, according to hypersonic theory including 
friction drag, can develop maximum lift-drag ratios of 4 and 
G (see Appendix B). We find (see Appendix C) that the 
L/D=A glide vehicle can radiate heat like a black body at a 
rate equal to the maximum average convective heat-transfer 
rate if the surface temperature is allowed to rise to about 
1 500° F. If the vehicle develops a lift-drag ratio of 6, then 
the allowable surface temperature must be increased to 
about 1800° F. These surface temperatures are high; 
nevertheless they are within the range of useful strengths of 
available alloys (see, e. g., ref. 7). Furthermore, they can, if 
necessary, be reduced somewhat by designing a less dense 
vehicle (or, more specifically, a vehicle of lower wing loading, 
W/S; see Appendix C). 

It is indicated, then, that the glide vehicle has the attrac- 
tive possibility of radiating back to the atmosphere a large 
fraction of the heat transferred to it by convection. As a 
result the mass of coolant required to protect the vehicle may 
be greatly reduced. Just as with the ballistic vehicle, how- 
ever (see ref. 3), it is evident that additional means, such as 
transpiration cooling, may be necessary to protect local hot 
spots on the surface, like the stagnation region of the nose. 
It is also well to note that the alleviating effects of radiative 
cooling are not limited to the glide vehicle alone, but would 
apply to any hypervelocity vehicle in level flight. 

We inquire now if the skip vehicle is capable of radiating 
heat at a rate comparable to the maximum convective heat- 
transfer rate. For this purpose it suffices to confine our at- 
tention to the first skip wherein the maximum convective 
heat-transfer rates are encountered (see eq. (64)). Calcula- 
tions of maximum average rates using equations (63) and 
(64), for long-range skip vehicles developing lift-drag ratios 
of 4 and 6, indicate that these rates are an order of magnitude 
higher than those for comparable glide vehicles. The cor- 
responding equilibrium surface temperatures of the skip 
vehicles are the order of two or more times as high as those of 
the glide vehicle. Accordingly they may be far in excess of 
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3000° F. We conclude, therefore, that long range, high L/D 
skip vehicles cannot radiate heat at a rate equivalent to the 
maximum convective rate because the surface temperatures 
required for radiation to offset convection would exceed the 
temperatures at which known structural materials retain 
appreciable strength (see ref. 7). 

Now the skip vehicle operating at lift-drag ratios in the 
neighborhood of 2 will absorb less heat than skip vehicles 
developing higher lift-drag ratios. However, as shown in 
Appendix C, the former vehicle still absorbs more heat than 
a comparable high -pressure-drag ballistic vehicle and it ac- 
crues no appreciable advantage by radiation. From the 
standpoint of heat transfer, then, it is indicated that the 
skip vehicle is inferior to both the ballistic and glide vehicles. 
That is to say, proportionately more coolant of one form or 
another would be required to protect the skip vehicle than 
would be required to protect ballistic or glide vehicles of the 
same range. The skip vehicle has other disadvantages as 
well. Certainly one of the most serious of these is the very 
high lateral loads (see fig. 3) that the vehicle would be re- 
quired to withstand during a skip from the earth’s atmos- 
phere. These loads, coupled with simultaneous high thermal 
stresses (due to high convective rates), would require the 
structure to be stronger and, consequently, heavier than that 
of a comparable glide vehicle. 3 For these and other reasons 
concerned with problems of stability, control, and guidance, 
the skip vehicle is thought to be the least promising of the 
three types of hypervelocity vehicle considered here. 

In essence, then, the preceding study has indicated that 
the ballistic vehicle exhibits the possibility of being relatively 
efficient for hypervelocity flight by virtue of the fact that 
aerodynamic heating can be markedly reduced through the 
artifice of using blunt, high -pressure-drag re-entry shapes. 
The disadvantage of using the relatively inefficient ballistic 
trajectory is counterbalanced by this advantage which tends 
to keep initial mass down by reducing coolant mass. The 
glide vehicle appears promising for hypervelocity flight 
because it has, coupled with the relatively high efficiency of 
the glide trajectory, the possibility of radiating a large 
fraction of the heat absorbed by 'convection. 

Up to this point we have considered the performance 
efficiency of the several types of hypervelocity vehicle by 
comparison with each other. It is of interest now to compare, 
insofar as is possible, the efficiency of flight of these vehicles 
with that of lower speed, more conventional type aircraft. 

COMPARISON OF HYPERVELOCITY VEHICLES WITH THE SUPERSONIC 

AIRPLANE 

In the analysis of powered flight it was found that the basic 
performance equation for hypervelocity vehicles could be 
written in a form analogous to the Breguet range equation. 
Thus, according to equations (7) and (9), we have for both 
hypervelocity and lower speed vehicles that 

fi=(s)/.Kln(i;) (84) 

3 Added weight means, of course, added coolant (see, again, eq. (61)) and one can easily 
demonstrate that ultimately the coolant is being added to cool coolant. This situation must 
obviously be avoided. 


where it is understood that the effective quantities are the 
same as the actual quantities in the case of the lower speed, 
more conventional aircraft. Now let us consider the product 
( L/D) e I e V e . Taking first the supersonic airplane we assume 
flight at a maximum lift-drag ratio of 6. The product I e V e 
for a ram-jet or turbojet can reasonably be expected to have 
a value of about 4.4X10 6 feet. 4 The product (L/D) e I e V e is 
then 26.4X10 6 feet for the airplane. Now let us compare 
these quantities with the corresponding quantities for a 
ballistic vehicle and let us presume that the range will be half 
the circumference of the earth. In this event, the effective 
lift-drag ratio for the ballistic vehicle is 2i r (see fig. 2) which is 
slightly greater than that for the airplane, while the effective 
velocity is just half the satellite velocity, or 13,000 feet per 
second. Let us again assume that the effective specific 
impulse is 300 seconds. In this case, the product of I e V e is 
3.9X10 6 feet and the product (L/D) e I e V e is about 24.5X10 6 
feet which is only slightly less than that for the supersonic 
airplane. Thus we have our first suggestion that the hyper- 
velocity vehicle is not necessarily an inefficient type vehicle 
for long-range flight. 

In order to pursue this point further, a performance effi- 
ciency factor (see eq. (10)) defined as 



has been calculated for ballistic and glide vehicles for / e =300 
seconds, and ranges up to the circumference of the earth. 
The corresponding quantity E has been calculated for the 
supersonic airplane (U>Ve=4.4X10 6 feet) for several lift-drag 
ratios. The results of these calculations are presented in 
figure 10 and we observe, as our example calculation sug- 
gested, that both the ballistic and glide vehicles compare 
favorably with the supersonic airplane for ranges in the 
neighborhood of and greater than half the circumference of 
the earth. The glide vehicle is again superior to the ballistic 
vehicle at lift-drag ratios in excess of 2 and, as a result, it 
compares favorably with the airplane at shorter ranges than 
the ballistic vehicle. 

It should be kept in mind, of course, that m f may be 
substantially greater than m P , the mass of the pay load. 
This point is significant because it reminds us that milm P , 
and not m^m/, is considered the better measure of cost. 
Thus, noting that milm P =(milm f )(m f /rn P ) y and recognizing 
that m f jm P is probably lowest for the ballistic vehicle, we 
anticipate that the ballistic vehicle would appear to better 
advantage than shown in figure 10. 

CONCLUDING REMARKS AND SOME DESIGN CONSIDERA- 
TIONS FOR MANNED HYPERVELOCITY VEHICLES 

During the course of this study it has been indicated that 
ballistic and glide vehicles can be operated at hypervelocities 
with the reasonable assurance that problems of aerodynamic 


* This value should hold approximately for any air-breathing engine— note that the maxi- 
mum value of I,V» is simply the product of the thermal efficiency (taken as 0.3— see, e. g., ref. 
8) and the specific heat content of the fuel (taken as 14.6X 10* feet for gasoline-type fuels). 
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Figure 10. — -Variation of performance efficiency factor with total range 
for ballistic and glide vehicles and the supersonic airplane. 


heating can be largely alleviated by proper design. Skip 
vehicles appeared substantially less promising in this as well 
as other respects. It was further demonstrated that on the 
basis of equal ratios of initial to final mass, the long-range 
hypervelocity vehicle compares favorably with the super- 
sonic airplane. These considerations suggest that the 
ballistic and glide vehicles have, in addition to the advantages 
usually ascribed to great speed, the attractive possibility of 
providing relatively efficient long-range flight. 

In view of these findings, it seems appropriate as a final 
point to consider what appear to be favorable design features 
of manned hypervelocity vehicles. It is fair to assume that 
the glide vehicle has the man-carrying capability if suitable 
living quarters are provided inside the vehicle, particularly 
as regards composition and temperature of the interior 
atmosphere. However, whether or not the ballistic vehicle 
has this capability is not obvious and requires some clarifica- 
tion. The principle question in this regard is the magnitude 
of the decelerations experienced by the vehicle and its 
occupants during atmospheric entry. Some light is shed on 
this matter by figure 11 where the maximum deceleration in 
g’s of a ballistic vehicle is shown as a function of range. 



Figure 11. — Maximum deceleration of ballistic vehicles during atmos- 
pheric entry. 


These decelerations were calculated bv the method of 
reference 3 for large entry angles, and with equations (13), 
(C L =0) for entry angles near zero, using the velocities and 
entry angles as a function of range given by equations (12). 
It seems reasonable to conclude from the results shown in 
figure 1 1 that the decelerations are in excess of those humanly 
tolerable except for very short range flight and for very long 
range flight. The latter case is of principal interest to us, 
and it is noted specifically that maximum decelerations can 
probably be kept to the order of 10 g ’ s or slightly less for 
ranges of the order of semiglobal and greater. It may be 
remarked further that decelerations exceed 5r/s for less than 
a minute, and they exceed 1 g for not more than about 3 
minutes. In this respect, then, (see ref. 9) the ballistic 
vehicle appears to be a practical man-carrying machine, 
provided extreme care is exercised in supporting the man 
during atmospheric entry. From the aerodynamic heating 
point of view the ballistic vehicle can, of course, be made 
especially attractive by employing the blunt body concepts 
of reference 3. 

The glide vehicle experiences maximum decelerations in 
g ’ s equal to approximately D/L (see eq. (31) and note 
D/W-+D/L as I T2 -*0), and so with any significant lift-drag 
ratio it is far superior to the ballistic vehicle in this respect. 
In addition, the glider has the important advantage of 
maneuverability during atmospheric entry. These factors 
and its potential for relatively high performance efficiency 
make the glider generally attractive as a man-carrying 
machine. 

It will be assumed that if the glider is to develop reasonably 
high lift-drag ratios it should be slender in shape. But the 
nose of the body and the leading edges of the wing (and tail 
surfaces) should be blunt to alleviate the local heating prob- 
lem. Blunting the nose of the body may not, if properly 
done, increase the drag of the vehicle (see refs. 10 and 11). 
Blunting the leading edge of the wing will, however, incur a 
drag penalty and thereby reduce the lift-drag ratio. This 
difficulty may be largely circumvented by sweeping the lead- 
ing edge of the wing. The contribution to total drag of the 
drag at the leading edge is, according to Newtonian theory, 
reduced in this manner by the square of the cosine of the 
angle of sweep for constant span. The question which arises 
is how does sweep influence heat-transfer rate. The nature 
of this influence (ref. 6) is shown in figure 12 and it is ob- 
served that sweep decreases heat-transfer rate very substan- 
tially, although not to the extent that it decreases drag. 
We are led then to the conclusion that the wing on a hyper- 
velocity glide vehicle which develops reasonably high lift- 
drag ratio should have highly swept leading edges. This 
observation coupled with the fact that wing weight should 
be minimized suggests for our consideration the low-aspect- 
ratio delta wing. In addition to the wing it is anticipated 
that a vertical tail will be needed to provide directional 
stability and control, and so we are led to imagine as one 
possibility a hypervelocity glider of the type shown in 
figure 13. 

The potential of the glider to have relatively high per- 
formance efficiency hinges strongly on the finding that the 
large majority of the heat convected to it may be radiated 
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Figure 12. — Effect of sweep on drag and heat transfer to circular 

cylinders. 



Figure 13. — Example high lift-drag ratio glider. 

away at reasonably low surface temperatures. But it is 
never possible to build a perfect radiation shield. There is 
always a certain amount of heat which leaks through the 
shield to the internal structure. As the duration of flight 
increases this heat leakage problem may assume major pro- 
portions if substantially more structure (or coolant) is re- 
quired to absorb the heat. If, at the same time, the action 
of aerodynamic forces has, at best, a minor influence on 
range then the high lift-drag-ratio glider may cease to be an 
attractive machine. For flights approaching global range 
these two factors tend to come into play. That is, flight 
time becomes relatively long (of the order of an hour and 
a half or more) with the attendant increase in seriousness of 
the heat leakage problem, while lift-drag ratio assumes a 



Figure 14. — Example high lift glider. 


relatively minor role in terms of performance efficiency (see 
fig. 10). Accordingly, it may be attractive to launch a global 
glider into a low altitude satellite orbit which it follows over 
the large majority of its range and from which it enters the 
atmosphere in the terminal phase of flight to glide the short 
remaining distance to its landing point. Under these cir- 
cumstances, the vehicle may be designed to minimize aero- 
dynamic heating during atmospheric entry and for this pur- 
pose we are attracted to the use of high lift 5 as well as low 
wing loading (see eqs. (76) and (77)) to reduce heating rates 
and surface temperatures. Accordingly, the vehicle may 
glide into the atmosphere at a high angle of attack for high 
lift coefficient, maintaining this attitude until speed has been 
reduced to a supersonic value where heating has become a 
relatively minor problem. The angle of attack may then be 
reduced to increase LID. thereby extending the glide and 
increasing maneuverability to achieve the desired landing 
point. For this type of application the vehicle might have 
more of the appearance shown in figure 14, again being of 
the delta-wing plan form but having a more or less rounded 
bottom and sides to minimize heating rates over the leading 
edge as well as the entire lower surface during re-entry. 
Such a configuration bears a resemblance to a motorboat 
and it may in fact be suited for landing on water as shown. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Dec. 10 , 1954 


5 High lift tends, of course, to mean increased decelerations because of reduced LID during 
atmospheric entry; however, even for L/D's of the order of unity these decelerations remain 
modest and they should not, therefore, constitute a serious piloting problem. 


APPENDIX A 


SIMPLIFYING ASSUMPTIONS IN THE 

' M' 


The assumption of small deflection angle (0<<1) was 
used throughout the study of the glide trajectory. In 
addition, equation (28) was simplified on the assumptions 
that 

(L/D)6« 1 (Al) 

and 


do L/DdV* 
2 ds 


(A2) 


The extent to which these assumptions are permissible can 
be checked by deriving an expression for ( L/D)0 and exam- 
ining its variation over a range of trajectory parameters. 

From equations (15), (30), and (31) the altitude of any 
point in a glide trajectory is found to be 



By retaining the assumption of small inclination angle, 
whereby 0~—dy/ds, and recalling that J= (2s/r 0 )/(L/D), we 
find the inclination angle by differentiating equation (A3). 
Performing this operation and making use of equation (30) 
reduces the expression for ( L/D)0 to 


L 0== _2_ r 1 _ 1 2.105X1Q-" 3 

D fir 0 1_] — (j — Ty)c y J V 2 


(A4) 


Since V 2 becomes very small toward the end of the trajectory, 
it is apparent from equation (A4) that the assumption of 
small ( L/D)0 cannot be justified in this portion of flight. 
The problem then is to determine the conditions under 
which (L/D)0 remains negligibly small over the major part 
of the trajectoiy. 

With the aid of equations (30) and (33), equation (A4) 
can be modified to the following form 


L/D 

$ 2 $ 




In * 


1 



(A5) 


For given values of L/D and total range equation (A5) 
determines the fractional part of the total range which cor- 
responds to a given value of (L/D)0. Since the deflection 
angle is always increasing, we can therefore determine the 
portion of the total range through which ( L/D)0 remains 
equal to or less than a given value. A computation of this 
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nature was performed for a value of {L/D)0< 0.05, and the 
results are presented in figure 15. From this figure we can 
see that except for short ranges and large lift-drag ratios, 
( L/D)0 (as well as 0) remains at a value less than 0.05 fojr 
better than 90 percent of the total range. 

The second assumption, equation (A2), can also be verified 
from the results of the analysis. By differentiation of equa- 
tion (A-4) we find that 


V 2 4F, 2 f (\~V/)e J ~l 

ds fir 0 2 (L/D) 2 L V 2 J 


(A6) 


while differentiation of equation (30) yields 


1 (T/n) ^ V.'( l-V, 2 )e J 

2 K1 ' ,U) ds “ r 0 


(A7) 


Dividing equation (A6) by equation (A7), and making use 
of equation (A4), we find that 


V2 d0 

ds _ 20 _ 2 \(r/n\fi i 

i (t it* [{ 1 m 


(A8) 


By comparing equation (A8) with the previous results ob- 
tained for ( L/D)0 (fig. 15), we can readily see that the assump- 
tion of equation (A2) is actually less stringent than that of 
equation (Al ) for values of L/D of the order of 1 and greater. 



Figure 15. — Portion of range where {LfD)e< 0.05 as a function of 
range for various values of lift-drag ratio. 
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THE RELATION BETWEEN 


CCS 

CdA 


AND 


The lift and drag coefficients for slender cones at small 
angles of attack can he expressed in the following manner: 

(\=2a (Bl) 

f d = ( ' d 0 A ot( \ (B2) 

from which it follows that the lift-drag ratio is 




a 


r, 


D r D +a(\ On+CJ/2 


(B3) 


It can he shown that equation (B3) has a maximum value 
when 

C d =CYI2 (B4) 

whereby 

{ < (, D)(L\D) ma =2C Di> (B5) 

Using equations (Bl) through (B5), one can express the 
maximum value of the lift-drag ratio in the following ways: 

(B6) 

The drag coefficient at zero angle of attack appearing in 
equation (B2) can be broken down into its component parts 
to yield 

( \=( C °X+(\S/A (B7) 

where (C D J p is the zero-lift pressure drag coefficient and 
C Fo is the zero-lift skin-friction coefficient based on wetted 
area. The skin-friction coefficient C Po in equation (B7) can 
be related to the equivalent skin-friction coefficient fV (see 
eq. (43)) by considering average conditions over the surface 
of the cone. Equating the friction drags as determined 
from free-stream and local average conditions, it is found 
that 


/~r //*y \ (pi)at>(Uf) c 

Lf o V F l)av p y2 


(B8) 


By referring to local average conditions on the body sur- 
face, the expression for C F ' , equation (40), can be written as 

r '—in \ av (V,) 

t F ~{C Fl ) ac ^ — (B9) 

Comparing equations (B8) and (B9) it is apparent that 


C F — Ci 


, (F,) 0 


(BIO) 


For slender shapes at hypersonic speeds, the local velocity 
does not differ appreciably from the free-stream value. 
Also, for small angles of attack, the skin-friction coefficient 


(I). 


FOR CONICAL MISSILES 


should remain fairly constant. Consequently, equation 
(BlO) can be written as 


('f= C f =C f ’ = const an t 
and equation (B7) then becomes 

/** ( C \ I C/S 


(B11) 


(B12) 


From equations (B5) and (Bl2) it can then be shown that 

['-¥-] « 

From the Newtonian impact theory, the zero-lift pressure 
drag coefficient for slender cones at hypersonic speeds can 
be expressed as 

(C d X=2S 2 (Bl4) 

where 8 is the semivertex angle of the cone. By further 
noting that for slender cones 


S 1 
A 5 


equation (B 12) then becomes 


C/ 


C d =28>+^ 

• 8 


(B15) 


(B16) 


For a given value of C/ it can be shown that equation 
(B 1 6) has minimum value when 6 


-~-Q¥T 

whereby, at 8 = 8 opt 


(B17) 


(B18) 


Obviously, then, the highest value of maximum lift-drag 
ratio (eq. (B6)) will be attained by the cone with the semi- 
vertex angle given by equation (B17). By substitution 
from equation (Bl8) into equation (B6), the optimum value 
of maximum lift-drag ratio is found to be 


(l\ — 1 ( 4 y*_ i 

\DJmax ^V3 VCf* / 2-y/3 8 0 


(BIO) 


By furtlier substituting the expression for minimum 
zero-lift drag coefficient, equation (Bl8), into equation 


# The remaining analysis assumes C F ' constant. Although this Is certainly not the 
situation in practice, the analysis provides an “ order of magnitude” estimate of pertinent 
parameters. 
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(B13), the following relation, corresponding to the condi- 
tion of optimum maximum lift-drag ratio, is obtained: 



(LID) 



(B20) 


With the aid of equations (B14) and (B16), equation 
(B13) can also be expressed in the following form, cor- 
responding to any maximum lift-drag ratio including the 
optimum value: 

From equations (B20) and (B21) it can readily be seen that 
in the case of the optimum (L/D) max 

^3=2 (B22) 

40 O pt 

from which it follows directly that 

f)’ <B23) 

so that equation (B21) may be written as 


(C/S\ 

| — 

1 

KcbA) 

(L/D) maz 2-|~| 

f-Y 

\&optJ 


With the aid of equations (Bl6) and (B23), the expression 
for any (L/D) max (eq. (B6)) can be shown to be 



(B25) 


and it follows directly from equation (B19) that the ratio 
rj of {L/D) max for any cone to that for the optimum cone is 


V= 


fippt 

8 



(B26) 


where rj is defined as the ‘‘lift-drag efficiency factor / 1 By 
substitution from equation (B24), the ratio of ( L/D) max to the 
optimum value can then be expressed in terms of ( C/S/C D A ) 
as follows: 


V= 



c/sy 

CoA) 


(B27) 


The dependence of rj on C F 'S/C D A is shown in figure 16. 
It should he noted, however, that for small values of rj the 
assumption of slender cones will be violated, although the 
results as shown will be qualitatively correct in that C F 'S/C D A 
will become exceedingly small for low values of ( L/D ) max , 
regardless of body shape. 



Figure 16. — Variation of lift-drag efficiency factor with drag param- 
eter for cones. 


APPENDIX C 


COMPUTATION OF HEATING ASSOCIATED WITH ROCKET VEHICLES 


RADIATION OF HEAT FROM GLIDE VEHICLES 

From equation (72), the maximum time rate of total heat 
input to the glide vehicle can be expressed as 

fdQ\ _ o (dH a \ mgVs ( C F 'S \ , pn 

\<U /max \ dt J mai 3-y/3(L/D) \CdA J 


The rate of heat radiation from the vehicle is given by the 
relation 


S 


dHs 

dt 


kT K *S 


(C2) 


Using equations (Cl) and (C2), the requirement for con- 
tinuous radiation of all convective heat input to a surface 
at a temperature 2000° R can be expressed as 


mg C/S 


S {L/D) C d A ~ 


< 1.20 


(C3) 


If a value of {L/D ) max — 6 is assumed, values of the para- 

Opt 

meter C/S/C D A and cone angle, 5, can be determined as a 


function of (L/D) max from the analysis given in Appendix B. 
A vehicle weight of 5000 pounds with a maximum diameter 
of 3 feet is assumed whereby equation (C3) can be evaluated 
for various {L/D) giving the results in the following table: 


LID 

Cr'S/CoA 

i, 

deg 

S, 

sq ft 

mg Cr*S 

S(L/D) CdA 

6 

0.333 

2. 75 

147 

1. 89 

4 

.0600 

6. 73 

60.2 

1.24 

2 

.00710 

14.3 

28.6 

.619 

1 

.000890 

29.6 

14.3 

.312 


.000115 

55.8 

8. 55 

.135 


We see, therefore, that at surface temperatures of 2000° R 
and for an L/D of 4, this glide-type vehicle can radiate heat 
at a rate equal to or greater than the maximum convective 
heat rate. 

RADIATIVE AND CONVECTIVE HEAT TRANSFER ASSOCIATED WITH SKIP 

VEHICLES 

In this section the problem is to determine the extent to 
which heat absorbed by a skip vehicle in the first skipping 
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phase can be reradiated during the subsequent ballistic 
phase. The quantity of heat absorbed in the first skipping 
phase has already been obtained in the heating analysis, 
(eq. (59) for n=l) 


Q ! 1 C/S/ 

\mVf 2 C d A \l-e L,D i 


(C4) 


where the total heat absorbed throughout the entire trajec- 
tory is 


Q t _ 1 C/S 
\mV f 2 2 CuA 


(C5) 


In order to determine the heat radiated, three quantities 
must be determined: 

1. Temperature of the vehicle at the start of the second 

ballistic phase 

2. Temperature of the vehicle at the end of the second 

ballistic phase 

3. The time duration of the second ballistic phase 

To determine the first of the above quantities, we employ 
the relation for heat absorbed 


Q x =cWAT (C6) 

where c is the specific heat of the material, W e is the effective 
weight of material absorbing heat, and AT 7 is the tempera- 
ture rise during the first skip. If it is assumed that 1/3 of 
the missile weight will absorb heat, equation (C6) becomes 



where m is the total mass of the vehicle. It is assumed that 
the material has a specific heat of 0.11 Btu/lb °R. If it is 
also assumed that the temperature at the start of the first 
skip is 500° R, equation (C7) becomes 

2Vr 3 =500+l.lX10 -3 (~~) (C8) 

which defines the temperature at the beginning of the 
second ballastic phase. 

To find the temperature at the end of the second ballistic 
phase, we equate the radiant heat-transfer rate from the 
body to the rate of heat loss in terms of the temperature 
drop of the body 

—JcT 4 S dt=cW e dT (C9) 

This expression can be integrated to yield 

TJ= r (CIO) 

(7.95X10" 15 )St+^ 1 

1 ex 2 

for a vehicle weight of 5000 pounds (effective absorbing 
weight of 1667 pounds) where T en<1 is the temperature at 
the end of the second ballistic phase and t is the total flight 
time of the second ballistic phase. The total heat lost by 
radiation can now be expressed in terms of the temperature 
drop as 

QR 2 ={T ex -T tn2 )W t C 

or 

Q R =lAl X 10 \T ex - T 9n2 ) (Cl 1) 


The time of flight in any ballistic trajectory can be shown to 
be 


,2 r 0 

l ~v s 


—ff/' [ ta ° *' + vI =2 tan " (Vra tan f)] 


where 




tan 0 


/ 


. $ $ 
sin — + tan 0 f cos ~ 


f =tan- 


fein Of cos 0\ 

= — cos 2 0 f i 

<V 2 f 


(Cl 2) 


The foregoing relations were applied to a computation of 
the radiative cooling of a missile weighing 5000 pounds and 
traversing a total range of 3440 nautical miles (<l>=1.0). 
Values of 6 f were obtained in the motion analysis, and values 
of C/SIC d A and S obtained in the previous calculation 
with regard to the glide missile will apply to this case also. 
The computations are summarized in the following table. 
Note that the case of L/D — }i is essentially the ballistic 
vehicle (see fig. 9) . 


LID 

er, 

deg 

V/ 3 

QiXIO-s 

Qi/Qr 

t, 

sec 

T«x j t 

°R 

T..„ 

°R 

Q*,X10-s 

Qk 2 /Qi 

6 

12.5 

0. 275 

3115 

0. 135 

213 

2710 

1490 

1725 

0. 554 

4 

17.0 

.315 

1470 

.258 

335 

1542 

1323 

316 

.211 

2 

24.0 

.525 

549 

.575 

395 

889 

885 

27 

.049 

1 

27.5 

.620 

122 

.853 

247 

587 

585 

3 

.023 


30.0 

.650 

19 

.985 

80.4 

514 

514 

0 

0 


We see, therefore, that the quantity of heat which must 
be absorbed by this skip vehicle decreases rapidly with 
decreasing lift-drag ratio. The quantity of heat which must 
be absorbed by a ballistic vehicle (L/Z)~ 1/2) is almost 
negligible compared with the quantities associated with 
vehicles with an L/D= 2 or greater. Comparison of the heat 
absorbed in the first skipping phase with the heat radiated 
in the second ballistic phase indicates no appreciable ad- 
vantage is obtained due to radiation for values of LID =2 
and lower. To be sure, this situation could be substantially 
altered (near LjD=2) by allowing the surface temperatures 
to reach higher values during the skip; however, it seems 
unlikely that the net heat absorbed by the skip vehicle 
could ever be reduced to the low value of the ballistic vehicle 
for any reasonable surface temperature. 
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